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Abstract — In this paper we consider prioritized maximal 
scheduling in multi-hop wireless networlis, where the scheduler 
chooses a maximal independent set greedily according to a 
sequence specified by certain priorities. We show that if the 
probability distributions of the priorities are properly chosen, 
we can achieve the optimal (maximum) stability region using 
an i.i.d random priority assignment process, for any set of 
arrival processes that satisfy Law of Large Numbers. The pre- 
computation of the priorities is, in general, NP-hard, but there 
exists polynomial time approximation scheme (PTAS) to achieve 
any fraction of the optimal stability region. We next focus on 
the simple case of static priority and specify a greedy priority 
assignment algorithm, which can achieve the same fraction of 
the optimal stability region as the state of art result for Longest 
Queue First (LQF) schedulers. We also show that this algorithm 
can be easily adapted to satisfy delay constraints in the large 
deviations regime, and therefore, supports Quality of Service 
(QoS) for each link. 

Keywords: Maximal scheduling, wireless networks, stabi- 
hty, delay, priority, large deviations. 

I. Introduction 

Efficient scheduling algorithms in wireless networks have 
been the subject of intensive research in the past few years. 
A fundamental issue is that the optimal scheduling algorithm, 
which can achieve a superset stability region of any schedu- 
ler, involves solving a Maximum Weighted Independent Set 
(MWIS) problem [1], which is, in general, NP-Hard [2]. This 
phenomenon is exaggerated by the requirement that the MWIS 
problem has to be solved in every time slot, which renders 
it impossible to be implemented in applications due to the 
high computing resource consumptions. As an alternative, low 
complexity suboptimal schedulers with guaranteed efficiency 
ratio are needed, among which, maximal scheduling has been 
the focus of recent research. Similar to the maximal matching 
in the switch scheduling literature [10], maximal scheduling 
is a low complexity algorithm in wireless networks with 
a constant efficiency ratio which is inversely proportional 
to the maximum number of independent links in a link's 
neighborhood [3], [5]. Not only is it sound theoretically, 
maximal scheduling is also well supported in practice by 
efficient distributed algorithms with constant overheads (e.g. 
see [4]). 

Furthermore, it has been observed that the performance gua- 
rantee for maximal schedulers can be improved significantly 



by considering specific maximal schedulers, since the class of 
maximal schedulers is quite broad. Recently, LQF scheduling 
has been shown to yield a much larger stability region than 
the worst case maximal scheduling. During scheduling, an 
LQF scheduler produces a maximal independent set greedily 
following a sequence according to the queue length order of 
each link, from the longest to the shortest. It is shown that in 
certain networks, where the topology satisfies the so called 
"local-pooling" condition [8], LQF is optimal. For general 
networks, the performance of LQF scheduling can be bounded 
by its "local-pooling factor" [9], which is a function of the 
network topology and interference model. In the worst case, 
LQF scheduling can guarantee an efficiency ratio between 1 /6 
and 1/3 in geometric networks with the K-hop interference 
model. 

Promising as it is, presently there are still some open issues 
for the LQF scheduling. First of all, the stability region of 
LQF is not well characterized, i.e., it is unclear whether a 
network is stable or not under an arrival process whose rate 
is outside the stability region of the the worst case maximal 
scheduler Moreover, unlike maximal matching for switch 
scheduling [10], where a centralized controller is available, 
distributed LQF scheduling in wireless networks is hard to 
implement, and is subject to performance degradation due to 
asynchronous queue length information updates. Finally, it is 
hard to adapt LQF scheduling to support QoS in wireless 
multi-hop networks, due to the intractability of its analysis, 
which allows its performance to be understood only in certain 
networks with sufficient symmetries [11]. 

Realizing the limits of LQF scheduling, we try to improve 
the performance of maximal scheduling from a different 
perspective. Note that in LQF scheduling, in each time slot n, 
the scheduler picks the links following a sequence according 
to queue lengths. We generalize this to prioritized maximal 
scheduling, where in each time slot the scheduler chooses a 
maximal schedule following a sequence specified by a priority 
vector p(n), which is not necessarily the same as the queue 
length orders. The priority process p(n) is, in general, a 
random process, and includes the priorities generated by the 
max-weight scheduling [1] and LQF scheduling as special 
cases. Therefore, it can achieve the optimal stability region, 
if the priorities are properly chosen. In this paper, we show 



that a simple i.i.d random process p{n) suffices to achieve 
the optimal stability region. Furthermore, the distribution of 
p(n) can be computed, or approximated by a PTAS, if the 
arrival rate value provided as of each link is available. In 
applications, the arrival rates can be obtained from either 
online estimation or the parameters by upper layer services, 
such as digital voice or encoded video. Since arrival rates 
parameters and pre-computation of the priorities are only 
needed when the network topology changes, this combined 
priority pre-computation and maximal scheduling approach is 
suitable for slowly changing wireless multi-hop networks. 

Next we consider a special class of prioritized maximal 
schedulers, where p(n) is constant over time. In principle, 
this is simpler than the LQF scheduler, where dynamic prio- 
rities are used, and its performance is not influenced by 
asynchronous queue length updates. Interestingly, we show 
that if the priorities are properly chosen, we can achieve the 
same efficiency ratio as LQF scheduling. Moreover, we give a 
specific lower bound characterization of the stability region 
of prioritized maximal scheduling with a fixed p. Further, 
we provide an algorithm that can successfully compute a 
stabilizing priority p as long as the arrival rate vector is 
inside the lower bound region of any fixed priority scheduler. 
Therefore, it is easy to check the stability of the prioritized 
maximal scheduler for any arrival process that is within a 
certain guaranteed fraction of the optimal stability region. 

Finally we consider the delay constrained scheduling pro- 
blem in wireless networks. With a major class of services 
being real time services, which are sensitive to congestion in 
the network causing buffer overflows, the issue of determining 
and guaranteeing delay in wireless networks is of equal, if 
not more, importance compared to the issue of stability. It is 
desirable to implement a scheduling algorithm such that the 
queue overflow probability for each link is below a certain 
small threshold. In this paper we try to solve this problem 
by using maximal scheduling with constant priority. We first 
analyze delay guarantees of the worst case maximal scheduler, 
as well as those of the prioritized maximal scheduler in the 
large deviations regime, and formulate an upper bound on 
the queue overflow probability. By exploring the similarity 
between the stability constrained and the delay constrained 
scheduling problems, we then propose a greedy priority 
assignment algorithm, which is adapted from the priority 
assignment algorithm designed for the stability case, with the 
guarantee that it will generate a satisfying priority p as long 
as the delay constraint can be achieved by some priority. 

The paper is organized as follows. Section HIl introduces the 
system model. Section |lll] describes the prioritized maximal 
scheduler with random priorities. We consider the constant 
priority case in Section IIVI and adapt it to support delay 
constraints in Section [V] Section [VTl concludes this paper. 

II. System Model 

A. Network Model 

We model the topology of the network as a directed graph 
G — (y, E), where V is the set of user nodes, and E is the set 
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Fig. 1 . (a) is an arbitrary network of 6 links, and (b) is its interference graph. 
The black nodes form a maximum independent set, while the gray nodes form 
a maximal independent set. 

of communication links. A link i = {u,v) ^ E only if node 
V is in the transmission range of node u. The interference is 
modeled by an undirected interference graph Gc — {VcE,,), 
where Vc is the set of links, and Ec is the set of pairwise 
conflicts. Two links {i,j) G Ec if and only if they are not 
allowed to transmit together For example, in the primary 
interference model, the only constraint is that a user node can 
not transmit and receive simultaneously. Therefore, two links 
E Ec if and only if they share a common node, in which 
case the scheduling problem is reduced to a matching problem 
in G. This model arises naturally in switch scheduling, and is 
suitable for wireless networks using Bluetooth or FH-CDMA 
physical layers [3]. If secondary interference is considered, 
one model is the K-hop interference model [2], which requires 
that two links within K hops can not transmit at the same 
time (Note that the 802.11 DCF (Distributed Coordination 
Function) corresponds to K — 2). See Fig.[T]for an illustration 
of the network topology and interference graph for a sample 
network consisting of 6 links. 

For each link i, define its neighbor set as Ni — {j : 
G Ec}- Thus, the transmission of link i is successful 
if and only if no link in Ni is transmitting. Denote the 
interference degree of link i as A^, which is the cardinality 
of the maximum independent set in the subgraph formed by 
{i} U Ni. It has been shown that A — max^gy^ A^ is related 
to the efficiency ratio of maximal scheduling [5]. To analyze 
the performance ratio of the prioritized maximal scheduling, 
we need the following definitions. Denote a node removal 
sequence as r = (ii, Z2, . . . , which is a permutation 

of (1. 2, . . . , and consider removing the nodes from Gc 
according to the sequence r. Thus, ik denotes the index of 
the fcth node that is removed. If a certain node i is removed 

(r) 

in some step k (i.e., ik = i), define 61 as the interference 

degree of link i when it gets removed in the subgraph at that 

step (i.e., without nodes zi,i2, . . . ,?fe-i). Define 

S ~ minmax(5, , 
ren ieVa 

where 11 is the set of nl removal (permutation) sequences. 
It has been shown that 1/6 is a lower bound on the local 
pooling factor for a network [9], and hence, the lower bound 



on the efficiency ratio of the LQF scheduling. We will show 
later that 5 is also a lower bound on the efficiency ratio of 
maximal scheduling with constant priority. 

B. Traffic Model 

We assume the network is time synchronized, and that the 
transmitter of each link i is associated with an external arrival 
process Ai{n), which is the cumulative number of packet 
arrivals during the first n time slots. We further assume that 
the arrivals happen at the end of each time slot, and in each 
time slot, the packet arrivals are bounded above by a constant. 
The arrival processes are also subject to the Strong Law of 
Large Numbers (SLLN), i.e., 

lim Ai{n)/n — a-i for alH G Vc, (1) 

n — *cxj 

with probability 1 (w.p.l), where the constant a; is the arrival 
rate for link i. Note that these assumptions on the arrival 
process are quite mild, as the arrivals can be dependent over 
time slots and also among different links. As a consequence, 
in this paper we only focus on single hop traffic, since it 
is straightforward to generalize these results to the multi- 
hop case (since departures from one queue may very well 
be arrivals to another queue, in our mild assumption). 

C. Scheduling Model 

In each time slot, a scheduler tt chooses an independent set 
for transmission. In this paper, we are interested in maximal 
scheduling, where a scheduled set of links has the property that 
no other link can be added without violating the interference 
constraints. Denote the family of maximal independent sets as 
Ai. We will also treat a maximal independent set m e 
as a column vector such that rn^ = 1 if i G m, and m; — 
otherwise, as long as there is no confusion. The queueing 
dynamics at each link can be expressed as 

Q,{n) = Q,{Q) + A,{n) ~ D,{n) (2) 

where Qi{n) is the queue length at link i at the end of time 
slot n and Di{n) is the cumulative departures during the first 
n time slots. Clearly, Di{n) must ensure that Qi{n) > 0. 

For prioritized maximal schedulers, during scheduling in ti- 
me slot n, an independent set is produced, as described above, 
guided by a priority vector p(n), where Pi{n) represents the 
priority of link i. We assume that link i has higher priority 
than link j if and only if pi (n) < pj (n) and that all priorities 
are distinct. Thus, the scheduler chooses departures following 
a sequence specified by p, from the highest to the lowest, and 
schedules link i for departure if it has a nonempty queue and 
no higher priority neighbor has already been scheduled. 

The throughput of a scheduler tt is represented by its 
stability region At^, which is the set of stable arrival rate 
vectors under tt. We define stability to be rate stability, i.e., 

lim Ai{n)/n~ lim Di{n)/n = at w.p.l (3) 

n — >oo n — 'oo 

for all i G Vc.lt has been shown that the max- weight scheduler 
in [1] can achieve the maximum stability region, ^max, which 
is the interior of the convex hull of the maximal independent 
sets in Gc (we do not consider the boundary points in this 
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paper). We begin our study of prioritized maximal schedulers 
by considering the stability region of prioritized maximal 
schedulers with random priorities, in the next section. 

III. Scheduling with Random Priorities 

In this section, we analyze the stability region of prioritized 
maximal schedulers when the priority vectors (p(rt)) form 
a random process. First of all, as a subclass of maximal 
schedulers, any prioritized maximal scheduler can support the 
following region [5]: 

Anun = {a : a,; + ^ aj < 1, 1 < i < \Vc\}. (4) 

In order to obtain a better bound, we need to explore the 
impact of the priority vectors. In fact, the stability region of a 
priority scheduler is closely related to the choice of priorities. 
For example, consider a star shaped interference graph, with 
center link 1 and n — 1 outer links. If we assign the priorities 
such that the center link 1 always has the lowest priority, then 
for any arrival rate ai + J2j=2 1' ™ arrival 

process which makes link 1 unstable. Therefore, the stability 
region coincides with Amin- As another example, consider the 
following choice of priorities: in each time slot, we first solve 
the MWIS problem 

m*(n) argmax^g^Q(n)^m. (5) 

Then, we assign all Unks with m* {n) = 1 the highest priorities 
(the order does not matter) while the links with m*{n) — 
are assigned the remaining (lowest) priorities (the order does 
not matter). Since m*(n) is an independent set, all links in the 
independent set have the local highest priority, i.e., the highest 
priority in that specific neighborhood. After scheduling, the 
scheduled links maximize the weighted sum (|5]l, and thus 
-4niax, in this example, can be achieved [1]. However, the 
computation of priorities is prohibitive, since (|5]l requires 
solving an NP-hard problem in every time slot. In fact, a 
much simpler prioritized scheduler can also achieve ^max, 
by utilizing only an i.i.d random process (p(ri)). 
We need the following lemma. 

Lemma 1: For any arrival rate vector a and an i.i.d random 
process (p(n)) of priority vectors, if there exists a S* = 
5*2, ... , Sn) where set Si C Ni, such that 

ai+ ^ < Pr(pi < pj,Vj G S'i), (6) 

a is stable under the maximal scheduler with process (p(n)). 

Proof: See in Appendix. ■ 

Thus, this lemma provide a lower bound on the stability 
region achievable by (p(n)). The following is the main 
theorem of this section. It uses Lemma 1 to show that any 
a G -4max can be stabilized by a properly chosen (p(n)). 

Theorem 1: In an arbitrary network with n nodes, for any 
a G y^max, there exists a stabilizing prioritized maximal sche- 
duler with an i.i.d process (p(n)). Furthermore, the support 
of (p(n)) consists of at most n + 1 elements. 

Proof: Since a is in an open set .4„iax, there exists an e > 
such that a+ee G Anax, where e = (1, 1, . . . , 1)-^. According 



to the Caratheodory theorem, a + ee can be represented as a 
convex combination of at most n+1 maximal independent 
sets, i.e.. 



n+l 



(fe) 



(7) 



k=l 



where 0^0 and e^9 = 1. We associate a priority vector 
pC^) with each m*^'^) such that if mf''' — I, i has the local 
highest priority, and other priorities are arbitrary. During the 
scheduling in time slot n, the scheduler chooses p*^'') with 
probability 9^- Therefore, we have 

Pr(p, <pj,Vj e TV,) > 
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and the stability follows from Lemma 1 by setting Si = iV^ 
for all 1 < I < n. ■ 

Intuitively, Theorem 1 shows that the traditional coloring ap- 
proach, which is used for scheduling in networks with constant 
traffics, is still applicable in the presence of stochastic packet 
dynamics. Furthermore, the hardness remains almost the same, 
since ^ is equivalent to a coloring problem. Note ^ is only 
executed during the pre-computation phase. Thus, compared 
to the max-weight scheduler, which essentially solves an NP- 
hard problem in every time slot, the burden on the prioritized 
maximal scheduler during the scheduling phase is significantly 
relieved. Finally, we show that there exists PTAS for the pre- 
computation of the priorities. 

Theorem 2: Given a network and interference model, if there 
exists a PTAS for the MWIS problem, then there exists a PTAS 
for computing the the stabilizing priorities. 

Proof: Suppose we have a + e'e e (1 — e)-4,nax for some 
e' > 0. The priorities can be obtained by solving the following 
optimization: 

min 



T 

e X 



subject to A/x ^ a + e'e, x ^ 



(8) 



where M is the matrix whose columns correspond to all m e 
M. This problem can be solved by binary search over t, where 
in each step we assume e^x = t and solve the following 

s 



max 



subject to 



(Mx), 



> s 



i<i< m 



e^x = i,xhO. (9) 
By strong duality, this can be solved by the dual problem 
mm /(A) 

subject to A ^ 0, e^A = 1 (10) 

where /(A) is the optimal value of the following problem 

(Afx), 



subject to 



E^^ 

e^x = t,x ^ 0. 



l<i< \VA 



It has been shown that an e-approximate solution for (|9]l can 
be obtained by solving 0(e) -approximate solutions for (fTTT i 
0(|Vc|(e^^ +log|Vc|)) times, see [16] for an approximation 
algorithm using logarithmic potential reduction. The only thing 
remains is to check that (fTTT i is a MWIS problem. Note 
that ( fTTT i is a linear program, with a simplex constraint, and 
therefore, the optimal value is attained at a vertex, i.e., 

/(A) = max w(A)^m (12) 

mGM 

where w^{\) = showing that /(A) finds a MWIS. Thus, 
the theorem holds. ■ 
For most network and interference models, there exists a 
PTAS for the MWIS problem, see [7] for an example using 
the graph partitioning technique in geometric graphs. Note that 
one can also achieve similar performance by using a PTAS 
for the max-weight scheduler during the scheduling [12]. 
Compared to their approach, our method has the advantage that 
the scheduling phase has low complexity, which is independent 
of approximation ratio, since the approximation algorithm is 
only executed during the priorities pre-computation phase. 

IV. Scheduling with Constant Priority 

In this section we consider a special class of prioritized 
schedulers, where p(n) = p is a constant vector, i.e., the 
priority for each link is fixed. This is certainly easier to 
implement, since it does not require a global coordination 
variable, specifying the priority vector applying in each time 
slot. We will focus on both of the priority pre-computation 
and the performance guarantees. 

A. Priority Assignment 

In order to search for the optimal priority, we need to 
analyze the associated stability region for an arbitrary priority 
vector p. By setting Si = {j € Ni : pi < pj} for each link 
I e Vc and apply Lemma 1, we get a lower bound 

fli + ajl{pi - pj) < 1 Vi G K (13) 

where l(-) is the indicator function 

1 if a; > 
otherwise 



l(x) = 



(14) 



Denote this stability region as Ap and the stability region 
achieved by any constant priority as ^ = ^pAp. We are 
interested in computing a stabilizing priority for a given 
a ^ A. We will first propose the algorithm, and then show 
its optimality. See also in [6]. Note that the proof below is 
new, compared to [6], and so is the condition that A^ A' . 



Algorithm 1 Stable-Priority (Gc,a) 



(11) 



for /c = 1 to \Vc\ do 

s ^ argminigy^{ai + Y^^eNi 
Ps n + I — k; 

Remove s from Vc and its incident edges from Ec', 
end for 
return p 



In each step, the algorithm chooses node s, the center node 
of a neighborhood, and assign it the next lowest priority. 

The following lemma is the key to the proof. 

Lemma 2: (Alternate definition of A) For the A defined 
above, we have A ~ A' where A is any set of rate vectors 
in M" that satisfies the following conditions: 

1) (Coordinate-convex) If a G and a' ^ a, then a' ^ A . 

2) (Worst-case stable) For any nonzero a G A, there exists 
1 < i < n such that > and + J2jeN- '^j ^ 

Proof: Suppose a G .4, then there exists p such that a G ^p. 
Thus, property (1) holds from the fact that is a poly tope, 
and thus contains any a' ^ a. For this priority p, denote i as 
the lowest priority link with positive arrival rate. Then, Si — (j>, 
so that ( fT3] ) is the same as property (2). Therefore we have 
AC A. 

Suppose a*^^' G A. From property (2) we have 

+ E 4' < 1 

for some ii. Assign pi-^ — \Vc\, and form a new rate vector 
a^^) from a^^^ by setting a-^^ = (which induces a reduced 
graph). From property (1), a*'^' is also in A, and from property 

(2) 

(2) there exists 12 ii such that aj^ > and 

We then set pi^ = |Vc| — 1 and repeat the similar procedure 
until the highest priority, i.e., priority (1) is assigned. Thus, a 
p vector is obtained. We claim that a G ^p, due to the fact 
that 

+ E ^AiP^-Pj) = Oif + E < 1 

for all 1 < k < \Vc\- Therefore we have A C A, and the 
lemma holds. ■ 

From the above proof, it is straightforward to prove the 
following theorem: 

Theorem 3: If a G ^, Stable-Priority will return a p such 
that a G ^p. 

B. Performance Guarantee 

We next analyze the performance of the maximal scheduling 
with the constant priority p generated by Stable-Priority. This 
can be measured by the efficiency ratio 7, which is defined as 

7 = SUp{cr : O-Anax G A}. 

To gain insight, we first consider the efficiency ratio of 
greedy graph vertex coloring, where a sequence of vertices 
are removed first and the colors are assigned in the reverse 
order. In fact, if the arrival processes have equal rate, and 
can be approximated by constant fluids, then the scheduling 
problem is reduced to a graph coloring problem, where the 
rate is inversely proportional to the number of colors used. 
Denote the removal sequence in a greedy coloring algorithm as 
^ — 42, • ■ • , in)- When the nodes are colored in the reverse 
order, we need at most rf'^'") + 1 = max^: d^^' + 1 colors, where 



I r ) 

dj. is the degree of node ik when it gets removed. Suppose 
that the above maximum is achieved at node ik*. When i^* 
gets removed, its neighborhood needs at least [ '^^(rj'"'" ] colors 
(in any optimal coloring), and at most d'^''^ + 1 colors in 
the sequential coloring. Therefore the efficiency ratio for this 
choice of sequence is lower bounded by l/6^'^\ In the worst 
case, we have the efficiency ratio of 1/A = 1/ maxrgn '^'■'"^ 
which corresponds to the efficiency ratio of the worst case 
maximal scheduler [5]. On the other hand, if we choose the 
sequence properly, we can achieve 1/6 — 1/ miuren <5''^-'- In 
the following we show that the efficiency ratio 1/S can, indeed, 
be achieved by A. 

Theorem 4: jAmax C A. 

Proof: Denote the removal sequence which achieves 5 as 
1" — «2, ■ ■ ■ , T'\Va\)- For any a G A define a sequence of 
arrival rate vectors as follows: 

a(i) = a 

aW = a('^-i)-a^^)e., 2<k<\VA 

where is the vector that is all-zero except for an 1 in the 
ifcth entry. In other words, a'*^') is obtained from a'^'''^^^ by 
setting the i^th entry to zero. We have for any k 

+ E «f ^ < ^ (15) 

due to the fact that at most 6 links's in i^'s neighborhood when 
it gets removed can transmit in each time slot. Now for the 
priority assignment such that pi^ = n + 1 — k for 1 < k < n, 
we have 

a., + E -Pj) = ^ + E 4'^ < ^ (16) 

for all ik. Thus we have j-a G Ap and the claim holds. ■ 
Interestingly, the efficiency ratio of 1/(5 has been shown to 
also be a lower bound on the local pooling factor [9] for an 
interference graph, which is the efficiency ratio promised by 
LQF scheduling. This is no coincidence, since in the special 
case of constant fluid arrivals with equal rate, both belong 
to the family of greedy coloring algorithms, and hence have 
similar performance guarantees. 

We next apply this result to certain networks and interfe- 
rence models to obtain examples of worst case guarantees. 
1 ) K-hop Interference: 

We first consider the X-hop interference model, which 
can be used for a large class of networks. For instance, the 
ubiquitous IEEE 802.11 DCF is usually modeled as a K- 
hop interference model with K — 2, due to the RTS-CTS 
message exchanges. In the i^T-hop interference model, two 
links G Ec if and only if the distance between one 

node in link i (transmitter or receiver) and one node in link 
j (transmitter of receiver) is less than a threshold Kr, where 
r is the transmission range of a node. Summarizing the result 
in [9], we have the following corollary. 

Corollary 1: In a geometric graph with A'-hop interference 
model, the prioritized maximal scheduling with constant prio- 
rity can achieve an efficiency ratio between 1/6 and 1/3. 



2) PHY-Graph: 

We next consider the PHY-Graph [17], which is a more 
reahstic interference model explicitly incorporating the physi- 
cal layer parameters, i.e., the signal-to-interference-plus-noise 
ratio (SINR). In the PHY-Graph model, two Hnks e Ec 
if and only if either of the following is true: 1) the distance 
between the transmitter of link i and the receiver of link j 
is less than clj, or 2) the distance between the transmitter 
of link j and the receiver of link i is less than ck, where 
c = {SNRt)i is a function of the SINR threshold SNRt and 
path loss exponent k, and li,lj are the link lengths of i,j, 
respectively. From their coloring bound, which is essentially 
an upper bound on S, we have the following corollary. 

Corollary 2: In a geometric graph with PHY-Graph inter- 
ference model, for fixed k, the efficiency ratio of prioritized 
maximal scheduling with constant priority is a nondecreasing 
function of SNRt- Particularly, when the SNRt is sufficiently 
high, the efficiency ratio is bounded above by 1/7. 

V. Delay-Aware Scheduling 

In this section we try to adapt the prioritized maximal 
scheduling to support QoS constraints, in the form of an upper 
bound on queue-overflow probability. Specifically, we assume 
that the system constraint on the queue over flow probability 
for link i is 

Pr(Qi(0) >B,)<e 

where Qi{0) is the stationary queue length and Bi is buffer 
capacity. For large Bi, we can assume that the buffer capacity 
if infinite and approximate the queue overflow problem as the 
following: calculate 9*, where 

0* = liminf logPr(Q,(0) > B,) > e'. (17) 

Here e = exp(— Sic') while 6* represents the delay exponent 
of Qi in the large deviations regime. In order to measure the 
delay performance of a scheduler vr in the large deviations 
regime, using a similar notation to the stability region Air, we 
define the delay region 

= {6* e M'; : 61 ^ 6*} 

where 6* is defined as Eqn. ( [TtI i. i.e., the set of guaranteed 
delay exponents under tt. In contrast to the stability region, 
where any the arrival process satisfying SLLN applies, the 
delay region is quite sensitive to the arrival process model, 
especially the burstiness of the process. Therefore, we have to 
first re-define the arrival process model before analyzing the 
delay region of maximal schedulers. 

A. Assumptions 

We need to slightly modify the arrival process model, so 
as to apply the Large Deviations Principle (LDP). We now 
assume that the arrival processes are independent among the 
links. For each link i, we will follow the model in [14]. 
Specifically, we have the following assumptions: 

Al: (The Gdtner-Ellis theorem applies) 



For each link i, the log moment generation function 

k^{e)= lim ilogE(e^^"(")) (18) 

exists for all 9, from which we can get the rate function via 
the Legendre transform 

A*(Ai) =sup(0A*-A,(0)). (19) 



A2: f77ie sample path LDP applies) 

For all s e N,ei,e2 > and for every scalars 
6o, 6i, . . . , hs-i, there exists iV > 0, such that for all n > N 
and all 1 ~ ko < ki <...< ks — n, 

> g-n'52+I]°Zo(fcj + i-fcj)A*(bj) ^20) 

where Pi is the probabihty of the following event: 

{\Ai{kj+i) - Ai{kj) - (fcj+i - kj)bj\ < ein,0 < j < s} 

Intuitively, this is the event that the arrival process is constrai- 
ned to lie within a tube around s linear segments of slopes 
bo, 6i, . . . , bs-i, respectively. 

A3: ( Convex dual analog of the sample path LDP) 

For all s G N, there exists > and a function g{-) with 

< 5(2^) < 00 when x > 0, such that for all n > A^ and all 

1 — ko < ki < . . . < kg ^ n, 

s 

E(exp(0^Z)) < exp{^(fc, - k,^,)A{9,) + g{9,)} 

where 9 = [9^,92, ... ,9,) and Z = (A,(fco), A,(fc2) - 

A, {ki), . . . , A,{ks) - Ai{k,^i)). 

It should be noted that the processes satisfying the above 
three properties form a broad class, which includes most 
common models for bursty traffic in realistic networks, such 
as renewal and Markov-modulated processes. For detailed 
discussions, see [13]. 

B. Worst Case Maximal Scheduling 

We first consider the delay region of an arbitrary maximal 
scheduler. In the following, we analyze the delay exponent for 
a fixed link i. Denote the original network with a chosen arbi- 
trary maximal scheduler as system S. Since it is analytically 
intractable to give an exact characterization of the delay region, 
we want to create a system S' such that 1) it is relatively easier 
to analyze the delay region in S' and that 2) the delay region 
of i in S' is a lower bound of that in S. Note that due to 
maximal scheduling, when link i has a nonempty queue, we 
can guarantee that, in each such time slot, there is at least one 
departure in i's neighborhood (which includes i). Therefore 
an obvious dominant system for 5 is a server with service 
capacity 1, which is fed by one queue of queue length equal 
to the sum of the queues in {1} U Ni. We can lower bound 
the delay exponent of Hnk i m S using the following bound 
based on S': 

Pr(g.(0) >x)< Pr(QKO) + ^ Q'j{0) > x). 



where the RHS corresponds to the single queue in S' (equal to 
the sum of queues in S) and can be calculated using standard 
large deviations techniques. However, empirical results show 
that this bound is quite loose, especially in the cases where 
some neighbors in Ni are likely to grow large, whereas Qi 
is not. To get a better lower bound, we need to analyze the 
queue length dynamics of link i alone. Since we are analyzing 
the worst case maximal scheduler, we assume a prioritized 
maximal scheduler where link i has the global lowest priority 
and the other links in Ni have the global highest priority, so 
that i can not transmit if any link in Ni has nonempty queue. 
Based on this, we create a system S' consisting of a clique 
formed by the nodes in {i} U Ni, and i is assigned the lowest 
priority. The priorities of the other links are arbitrary. We will 
show that S' is a dominant system over S in the following 
lemma. Intuitively, we are assuming the worst case correlations 
among the neighbors in Ni, such that the available time slot for 
Unk i is minimized. This, together with the scheduler which 
assigns link i the lowest priority, can be shown to dominate 
the original system S with any scheduler We will prove this 
domination in the following lemma. 

Lemma 3: For any a; > and n > 0, if Qi{Q) ~ Qi(0), we 
have 

Py{Q^{n)>x)<Pr{Q'^{n)>x) (21) 

where Qi{n), and Q'i{n) are the queue lengths of link i at 
time slot n ia S and S', respectively. 

Proof: We will show that in an arbitrary sample path u 
with the same arrival processes and initial queue lengths for 
both systems, we have Qi{n,Lu) < Q^{n,Lu) for all n > 0. 
We assume both systems use the first-in-first-out policy, and 
denote the arrival time and departure time of the fcth packet 
at link i in as tk and Sk, respectively. Similarly for S' we 
use t'f. and sj.. From the queueing equations in both systems, 
we have 



Q^{n) = Q»(0)+^l[o,„](ife)-^l[o,„](sfe) 

fc=l k=l 
oo oc 



k=l 



fc=l 



where Is (a;) is the indicator function that x G S*. Note that for 
any sample path oj, Qi{0,uj) — Qi{0,uj) and tk{uj) ~ t'f.{Lu) 
due to the assumptions. Thus to show that Qi{n) < Qi{n) 
w.p.l for all n > 0, it is sufficient to show that Sh{uj) < s'j,(ijj) 
for all fc > 1 and u. For simplicity we will drop the index u 
in the following and proof by contradiction. 

We will use induction on fc. Suppose that this is not true for 
the first packet, i.e., si > s[. In S', due to maximal scheduling, 
each time slot from t[ + l to — 1 is occupied by the neighbors 
of i (which is why the packet in i has not been served yet). 
Note that we always have 

D.,{n)+J2D,{n)>D[{n)+Y.^>^^ (22) 

for any n>l since both systems are work-conserving and the 
links in S' have edges between them (they form a clique), even 



when S does not. Thus from si > s'l, we have Di{s'i — 1) = 
D'^{s[ - 1) = and 

jeNi jeNi 
following Eqn. ( l22l l. Therefore, we have 

<"^^g;(s;-i) = J] 4^-1)- ^ i?;(,si-i) 

ieNi jeNi 

> 5: - 1) - V 1) 



= 



(fc) 

1) > 0. 



Equality (a) is explained as follows. Since i transmits in S' at 
time s\, no other neighbor is contending for that slot (neigh- 
bors have higher priority). It must be that all the neighbors 
have empty queues in the previous slot. The inequality (6) is 
explained as follows. Since i does not transmit in time slot 
s\ in S (we know it transmits at ,si > s'^, that must be 
due to some transmitting neighbor occupying that slot. Since 
arrivals occur at the end of each slot, that neighbor must have 
a nonempty queue in the previous slot s\ — \. Since this is a 
contradiction, we have proved the case for fc = 1. 

Now suppose it is true for up to fc — 1 packets. For the fcth 
packet, any time slot between + l to s'j. — 1 is occupied by the 
neighbors in S' . If Sk > sj,, we have Sk-i < s'j^^-^ < s'^. < s^, 
and 



(23) 



A(4 - 1) = - 1) 



1 



(24) 



Similar to (1231 ). we get 

0- E 



^ E 

= EQ^(4-i)>0' 

i.e., a contradiction. Thus the lemma holds. ■ 
Having shown the the dominance of S', we next analyze 

the delay performance of link i in S' in the following lemma. 
Lemma 4: For link i in system S' , we have 

lim — \ogVv{Q\{t) >x)^ e* (25) 

X — *oo X 

where 6* is the largest root of the following equation 

Aiie) + inf [ V A,(u) - m] = (26) 
o<u<e 



Proof: We construct a system S" equivalent to S', consisting 
of two nodes, where node 1 corresponds to link i and node 
2 corresponds to the links in A'^;. Thus, the queue of node 2 
is the sum of the queues of link i's neighbors. We assume 
that the arrival process to node 2 is equal to the sum arrivals 
in Ni, and that there is an edge between 1 and 2. During 
the scheduling, we assign higher priority to node 2. Note that 
from link i's perspective, S' and S" yield the same queueing 
dynamics. 



For a 2-node system, according to [14], the delay exponent 
0* for queue 1 is given by the largest root of the equation 

Ai(0)+ inf [A2(u)-u]^0. (27) 

0<u<9 

Thus the lemma holds from the fact that the arrival process of 
node 2 is A2{n) = J2jeNi ^ji"^)- ■ 
From the above discussions, we have the following theorem: 
Theorem 5: In a network where the arrival processes satisfy 
A1-A3, we have Qa\u\ Q ©tt for any maximal scheduler tt, 
where Q^ain = {0 £ M+ : 6* -< 6**} and 6* is the largest root 
of Eqn. dllll. 

This theorem provides a lower bound on the delay region 
of any maximal scheduler. 

C. Delay-Aware Maximal Scheduling 

Having obtained the asymptotic limit of the worst case 
maximal scheduler, we next consider improving the guaran- 
tees through proper priority assignment. In the following we 
assume that traffic parameter A; [9) for each link i is available. 
Ai(0) = fli, so that h.i{6) is more detailed than a^. First we 
generalize Theorem 5 to prioritized maximal schedulers. 

Lemma 5: The delay region 9p {0 £ R" : 6* ^ 0*} can 
be achieved by the prioritized maximal scheduler tt^, where 
9* is the largest root of the following equation 

k^{9) + inf [ y kj{u)l(j), -pj)-u]={) (28) 

Proof: Note that the neighbors in Ni with lower priorities 
than link i are invisible to link i. Therefore, we create a 
dominating system S' consisting of a clique formed by link 
i and its neighbors that have higher priorities than link i. 
Using a similar argument, one can show that Lemma 3 holds. 
Therefore, Lemma 5 holds, following Lemma 4. ■ 

Similar to A, we denote Q = Up0p, i.e., the delay 
region guaranteed by any constant priority scheduler We are 
interested in computing a proper priority p, when given a QoS 
constraint in the form of the vector 9. We will show this can 
be solved by an algorithm adapted from Stable-Priority, under 
the assumption that e 0, by showing that Q has a structure 
similar as A. 

Lemma 6: 8 = 0', where Q' is any delay region in M" 
that satisfies the following conditions: 

1) 0' is coordinate-convex; 

2) For any nonzero 9 G 0', there exists 1 < i < n such 
that 9i > Q and 9i is less than the largest root of the 
following equation 

K{9)+ inf [S^ KAu)l{9.)-u]^Q (29) 

j^N, 

where l(-) is the indicator function in (fT4l i. 

Proof: Suppose 6* G 0. Then, 9 G 0p for some p. Therefore 
property (1) holds following the fact that 0p is coordinate- 
convex. We can assume that p is such that the links with zero 
delay exponents have the lowest priorities. Next, consider the 



lowest priority link i with 0^ > according to p. We have 
9i < 9* where 9* is the largest root of the following equation 

U9.,) + inf [ y A,{u)l{p, ~Pj)-u] = 0. (30) 

Thus property (2) holds from the fact that l{pi —pj) = l(^^j) 
for all j G Ni. Hence we have C 0'. 

Now suppose G 0'. We first assign the lowest priorities 
to any link i such that 0^^^^ = 0. According to property (2), 
we can find link ii such that < 9- < 9) where 9) is 
described in Eqn. ( l29l l. We assign ii the current lowest priority 
available and set 9^"^^ = to get a new vector 9^'^^ G 0. Repeat 
the above process until all the links are assigned priorities. 
Denoting the resulting priority vector as p, it is easy to check 
that 6''^^^ G 0p, and hence 0' C 0. Therefore the lemma 
holds. ■ 

Based on the above lemma, we can construct a priority 
assignment algorithm, such that whenever 6* G 0, the al- 
gorithm will output a satisfying priority. This is the delay- 
equivalent of Stable-Priority algorithm, which was constructed 
using Lemma 2. 

Algorithm 2 Delay-Priority (G'c,6') 
1: for k ~ \ to n do 
2: for « G do 

3: Compute the largest root 9* of the equation: 

Ki{9) + info<„<eEj6Ar, A, (m)1(0j) - u] = 
4: end for 

5: s ^ min{z G K : 9, < 9*} 
6: Ps ^ 71 + 1 — fc 
7: 9s 

8: end for 
9: return p 

We conclude with the following theorem, which is the 
delay-analogue of Theorem 3. 

Theorem 6: If 6* G 0, Delay-Priority will generate a priority 
p such that the G 0p. 

VI. Conclusion 

This paper considered the prioritized maximal scheduling 
problem in multi-hop networks for arbitrary correlated arrival 
processes. We first considered the random priority case and 
showed that one can achieve the optimal stability region 
with i.i.d priorities. Then we focused on the constant priority 
case and proposed a priority assignment algorithm, which, 
combined with maximal scheduling, can achieve an efficiency 
ratio, which is the same as that achieved by the state of art 
LQF scheduling. We also analyzed the delay performance of 
maximal schedulers in the large deviations regime, assuming 
independent arrival processes, and proposed a delay-aware 
prioritized maximal scheduling algorithm. 

Appendix 

Proof of Lemma 1: Due to space limit, we only give a 
outline of the proof, for a detailed discussion about fluid limits 
please refer to [5] and [15]. 



1) The existence of fluid limits: Note that the support of 
functions Ai{n),Di{n) and Qi{n) is N, we extend it to 
R+ using Hnear interpolation, which results in continuous 
functions (e.g., Ai{t), t G K+). For any sample path w, define 
a family of functions as 



lim 

k — >oc 



(31) 



where /(•) could be Ai{-),Di{-) or Qi{-). Since both the 
arrivals and the departures in each time slot are bounded, and 
by linearity in (|2|i, the functions defined in dSTT l are Lipschitz 
continuous with the same Lipschitz constant (irrespective of 
r), and hence equi-continuous on any compact interval. Since 
they are also uniformly bounded on [0, t], according to Arzela- 
Ascoli theorem, in any [0,t], there exists a subsequence r„j. 
and continuous functions Ai{-), Di{-) and Qi{-), such that 
(uniform convergence) 

lim sup |/["'=(r,w) - /j(T,tj)| = 

fc^oo re[0,t] 

where _/(•) _ could be A{-),D{-) or Q{-). Any 
{Ai{-), Di{-),Qi{-)) satisfying the above is defined as a 
fluid limit. 

2) Properties of fluid limits: For any fluid limit, we have 
Ai{t) = Qit w.p.l, since Ai{n) satisfies SLLN, and 

Q^{t) = Q^{0) + a^t - D,{t) (32) 

for all i G Vc. The network is stable if, for any fluid Umit with 
0,(0) = 0,Vi, we have Q,{t) = for all i and all t > 0. 

3) Stability proof of Lemma 1: Define the Lyapunov function 

To prove the stability, it is sufficient to prove that if Li{0) ~ 
for all i, we have Li{t) = for all i and t > 0. Suppose that 
this is not true, then there exists i € Vc and t > such that 

L,{t) = max Z,(t) (34) 

Te[o,t] 



Q^{t) = .T > 0, 



(35) 



for details see the proof for (25) and (26) in [5]. Since Qi{t) 
is a uniformly continuous function, there is a t' < t such that 
Qi{t) > x/2 for T e [t',t]. Note that since Qi{t) is a fluid 
limit, there is a subsequence r„j. such that 



Qi[T) = hm > - 



(36) 



for every re [t' ,t]. Particularly, for large enough k, we have 

Q^{rn,T) > ^r„, > 1 Vr e [t',t] (37) 

i.e., the queue of link i is nonempty during [r„^t', r„j.t]. 
Therefore we have 

je{t}u{N,/Si) 
]e{i}uiNi/s,) 

Due to the randomized priority scheduling, since queue i is 
non-empty, if i has higher priority than any link in Si in a 
certain slot, we can guarantee one departure in {i} U {Ni/Si). 
Since the priorities are i.i.d across time slots, by SLLN 



> {t — t') Pr(queue i has local highest priority) 
= Pr(p, <p„Vj e 5,)(t-t'). 
Note that if, 

k 'OC> '^Uk 

= K+ E (SLLN) 

then we have 

U{t) - U{t') = lim 4"'' [t) - [t') < 

k — > oo 

which contradicts Eqn. ( [34l l. Therefore the stability of the 
network is proved and the lemma holds. ■ 
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